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£_| . Abstract 

O 1 

I We present a list of equivalent expressions and extensions of Bougerol's cele- 

brated identity in law, obtained by several authors. We recall well-known results 
and the latest progress of the research associated with this celebrated identity in 
many directions, we give some new results and possible extensions and we try to 
point out open questions. 



O 



CLh ■ AMS 2010 subject classification: Primary: 60J65, 60J60, 60-02, 60G07; 

i secondary: 60G15, 60J25, 60G46, 60E10, 60J55, 30C80, 44A10. 

Key words: Bougerol's identity, time-change, hyperbolic Brownian motion, subordina- 
tion, Gauss-Laplace transform, planar Brownian motion, Ornstein-Uhlenbeck processes, 
two-dimensional Bougerol's identity, local time, multi-dimensional Bougerol's identity, 
(*C) " Bougerol's diffusion, peacock, convex order, Bougerol's process. 

OO; 

Contents 

o: 

(NJ . 1 Introduction 

2 Extensions of Bougerol's identity to other processes 4 

2.1 Brownian motions with drifts I 

5h i 2.2 Hyperbolic Brownian motion 



Bougerol's identity and subordination 

3.1 General results 

3.2 Bougerol's identity in terms of planar Brownian motion 

3.3 The Ornstein-Uhlenbeck case 



*Laboratoire de Probability et Modeles Aleatoires (LPMA) CNRS : UMR7599, Universite Pierre et 
Marie Curie - Paris VI, Universite Paris-Diderot - Paris VII, 4 Place Jussieu, 75252 Paris Cedex 05, 
France. E-mail: stavros.vakeroudis@upmc.fr 

^Probability and Statistics Group, School of Mathematics, University of Manchester, Alan Turing 
Building, Oxford Road, Manchester MI 3 9PL, United Kingdom. 



1 



Multidimensional extensions of Bougerol's identity \w 

4.1 The law of the couple (sinh(/3 t ), sinh(L i )) 10 

4.2 Another two-dimensional extension 12 

4.3 A three-dimensional extension 14 



The diffusion version of Bougerol's identity 

5.1 Bougerol's diffusion 

5.2 Relations involving Jacobi processes 



16 
IS 



6 Bougerol's identity and peacocks [19 



7 Further extensions and open questions [20 



A Appendix: Tables of Bougerol's Identity and other equivalent expres- 
sions 



A.l Table 

A. 2 Table 

A. 3 Table 

A. 4 Table 

A.5 Table 



Bougerol's Identity in law and equivalent expressions (u > fixed) 
Bougerol's Identity for other 1-dimensional processes (u > fixed) 
Bougerol's Identity in terms of planar Brownian motion (u > fixed) 

Multi-dimensional extensions of Bougerol's Identity 

Diffusion version of Bougerol's Identity (relations involving the Ja- 



21 
21 
22 
23 



cobi process) [25 



Bibliography [25 



1 Introduction 

Bougerol's celebrated identity in law has been the subject of research for several authors 
since first formulated in 1983 [Bou83j. A reason for this study is on the one hand its 
interest from the mathematical point of view and on the other hand its numerous appli- 
cations, namely in Finance (pricing of Asian options etc.)-see e.g. |Yor92[ IDuf00[ lYorOlj . 
However, one still feels that some better understanding remains to be discovered. 

This paper is essentially an attempt to collect all the known results (up to now) and 
to give a (full) survey of the several different equivalent expressions and extensions (to 
other processes, multidimensional versions, etc.) in a concise way. We also provide a 
bibliography, as complete as possible. For the extended proofs we address the reader to 
the original articles. 

Bougerol's remarkable identity states that (see e.g. |Bou83[ IADY97] and |Yor01| 
(p. 200V), with (B u ,u > 0) and (f3 u ,u > 0) denoting two independent linear Brownian 
motional, we have: 

for fixed t, sinh^) (/ = ) f3 At {B) , (1) 

where A U (B) = J™ dsexp(2B s ) is independent of ((3 u ,u > 0). For a first approach of 
([T]), see e.g. the corresponding Chapters in |ReY99| and in |ChY12| . In what follows, 



§When we simply write: Brownian motion, we always mean real-valued Brownian motion starting 
from 0. For 2-dimensional Brownian motion we indicate planar or complex BM. 
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sometimes for simplicity we will use the notation A u instead of A u (-). 

Alili, Dufresne and Yor [ADY97J obtained the following simple proof of Bougerol's identity 

©: 

Proof. On the one hand, we define St = sinh(5 t ); then, applying Ito's formula we have: 

S t = [ VTTSI dB s + ~ f S s ds . (2) 
Jo * Jo 

On the other hand, a time- reversal argument for Brownian motion yields: for fixed t > 0, 

f>MB) = f z B *dls { = ] e B < f e~ B *d ls ee Q t , (3) 
Jo Jo 

where (7^, s > 0) denotes another 1-dimensional Brownian motion, independent from 
(B s ,s>0). 

Applying once more Ito's formula to Q t , we have: 

dQt = \Qtdt + (QtdB t + d lt ) = X -Q t dt + x/qT+1 d5 t , (4) 

where 5 is another 1-dimensional Brownian motion, depending on B and on 7. From (j2J) 
and 01]) we deduce that S and Q satisfy the same Stochastic Differential Equation with 
Lipschitz coefficients, hence, we obtain ([T|). □ 

With some elementary computations, from (JT]) (e.g. identifying the densities of both 
sides, for further details see [VakthlT] IBDY12a] ) . we may obtain the Gauss-Laplace trans- 
form of the clock A t : for every iGK, with a(x) = argsinh(a;) = log [x + \/T+ x 2 ) 



E 



1 ( x 2 M a'(x) ( a 2 (x)\ 

ex P KT- ■ ( 5 ) 



cxp 



At 1 \ 2A 



sft " V 2t 



where a'(x) = (1 + a; 2 ) . 

For further use, we note that Bougerol's identity may be equivalently stated as: 

sinh(|^|) {l = ] \(3\ Au(B) . (6) 

Using now the symmetry principle (see |And87] for the original note and |Gal08] for a 
detailed discussion): 

smh(B u ) l = (3a u (b), (7) 

where, e.g. B u = sup < s < u B s . 

In the remainder of this article we give several versions and generalizations of Bougerol's 
identity ([T]). In particular, in Section [2] we give extensions of this identity to other pro- 
cesses (i.e. Brownian motion with drift, hyperbolic Brownian motion, etc.). Section [3] is 
devoted to some results that we obtain from subordination and some applications to the 
study of Bougerol's identity in terms of planar Brownian motion and of complex-valued 
Ornstein-Uhlenbeck processes. In Section H] we give some 2 and 3 dimensional extensions 
of Bougerol's identity, first involving the local time at of the Brownian motion B, and 
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second by studying the joint law of 2 and 3 specific processes. In particular, in Subsection 
14.21 we give a new 2-dimensional extension. In Section [5] we generalize Bougerol's iden- 
tity for the case of diffusions, named "Bougerol's diffusions", followed by some studies 
in terms of Jacobi processes. Section deals with Bougerol's identity from the point 
of view of "peacocks" (see this Section for the precise definition, as introduced in e.g. 
[HP RYll] ). In Section [7] we propose some possible directions for further investigation of 
this "mysterious" identity in law with its versions and extensions and we give an as full 
as possible list of references (to the best of author's knowledge) up to now. Finally, in the 
Appendix, we present several tables of Bougerol's identity and all the equivalent forms 
and extensions that we present in this survey. These tables can be read independently 
from the rest of the text. 

We also note that (sometimes) the notation used from Section to Section may be 
independent. 

2 Extensions of Bougerol's identity to other processes 

2.1 Brownian motions with drifts 

Alili, Dufresne and Yor, in |ADY97 ] . showed the following result: 

Proposition 2.1. With fi, v two real numbers, for every x fixed, the Markov process: 



for every t > 0, has the same law as (sinh(ly ), t > 0), where (Y t ,t > 0) is a 
diffusion with infinitesimal generator: 




(exp(B t + fit)) ( x + / exp (-(B s + /is)) d((3 s + us) 



(8) 



n 




(9) 



starting from y = argsinh(a;). 

Proof. It suffices to apply Ito's formula to both processes and sinh(y^'^). 



□ 



It follows now: 



Corollary 2.2. For every t fixed, 




(10) 



In particular, in the case \i = 1 and v = 0: 





with e denoting a symmetric Bernoulli variable taking values in { — 1,1}. 



4 



Remark 2.3. With \i = — 1/2 and v = 0, we have that sinh (y£ 1 / 2 '°- ) j is a martingale. 



Indeed, with Y t = y} 1 ^ 2,0 ^ ; ltd 's formula yields: 



sinh(r t ) = / cosh(F s ) dY s + ]- [ sinh(F s ) ds 
Jo 2 Jo 

= f cosh(T s ) dB s - - tanh(y s ) ds +-/ sinh(F s ) ds 
Jo L ^ \ 2 J 



t 

cosh(Y s ) dB s . 

Jo 

Hence: 

M t = sinh(F t ) = /3 /o * ds ( cosh 2 (ys)) ee /3 /o ^ s ( 1+sinh 2 (ys)) , (12) 
and for this Markovian martingale, we have: 

M t = sinh(F 4 ) = f cosh(Y s )dB s = [ y/l + M 2 s dB s . (13) 
Jo Jo 

It can also be seen directly from (Q|) that {^X[~ l ^ 2,Q \t > 0^ is the product of two orthogonal 
martingales. This property is true because: 

B 



Y-(-l/2,0) _ f^i 

R u 



^ • (14 > 



with A[ = Jq ds exp(2i?i 1/ - ) ), (B\ u \t > 0) denoting a Brownian motion with drift, and 
(Rt,t > 0) a 2-dimensional Bessel process started at 0. Further details about this ratio 
are discussed in Sections [5] and We also remark that, with the notation of Section [ZJ 

2.2 Hyperbolic Brownian motion 

Alili and Gruet in |A1G97| generalized Bougerol's identity in terms of hyperbolic Brownian 
motion: 

Proposition 2.4. We use the notation introduced in the previous Subsection, that is: 
{Rt,t > 0) is a 2-dimensional Bessel process with R = and wee also denote by H an 
arcsine variable such that B^ u \ R and 5 are independent. Let <p be the function defined 
by: 

4>{x,z) = \j2e x cosh(,2) — e 2x — 1, for z > \x\. (15) 
Then, for fixed t, we have: 



/V> (/ = } (2E - 1)0 B?\ J R ? + (B^) 2 ) . (16) 



In particular, with v = 0, we recover Bougerol's identity: 

At {l = ] (2E - 1)0 (B U v^TT^) {l = ) sinh(^). (17) 
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This is an immediate consequence of the following: 
Lemma 2.5. (i) The law of the functional is characterized by: for all u > 0, 



E 



exp ( -— A\ 



vH l 2 \ dx e vx 



di -J== e~ z2t/2 J Q (u<P(x,z)), (18) 



<\x\ V2^t? 

where Jo stands for the Bessel function of the first kind with parameter \Leb7S^ . 
(ii) In particular, taking v = 0, for u > and x G K we have: 



x 



[--)E 



u 



exp ( — — A t }\B t = x 



+oo 



kl 



dz - e- z2t/2 J (u(f)(x,z)). (19) 



Proposition 12.41 follows now immediately from Lemma 12.51 by using the classical rep- 
resentation of the Bessel function of the first kind with parameter (see e.g. |Leb72] ): 



+i 



dr 



W-i vT 



cos(2x), 



(20) 



and remarking that (with H denoting again an arcsine variable), for all real £: 

J (0 = ^[exp(ie(2S-l))]. (21) 

Proof. ( Lemma \2.5\) 

With and denoting the modified Bessel functions of the first and the second kind 
respectively with parameter p = \J p 2 + v 2 (for p and v reals), we define the function 
: M 2 -> R + by: 



G^(u,v) 



2I^(u)K fl (v), u<v; 
2I fM {v)K fM {u), u>v. 



(22) 



First, using the skew product representation of planar Brownian motion, the following 
formula holds (for further details we address the interested reader to [A1 G97] ): 



(If oxi) j -— t ) E 



u 



exp -— A 



+oo 



dy e uy G^u,ue y ). (23) 



Using the integral representation (see e.g. |Leb72j . problem 8, p. 140): 

I^K^y) = - dr e'^Jo I ^2 cosh(r )xy — x 2 — y 2 \ , y > x. (24) 

1 Jlodv/x) v J 



'log(v/x) 

we can invert (123]) in order to obtain part %) of Lemma 12.51 
Part ii) follows with the help of Cameron-Martin relation. 



□ 
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3 Bougerol's identity and subordination 



In this Section, we consider (Z t = X t + iY t ,t > 0) a standard planar Brownian motion 
(BM) starting from x + iO, Xq > (for simplicity and without loss of generality, we 
suppose that x = 1). Then, a.s., (Z t ,t > 0) does not visit but it winds around it 
infinitely often, hence 9 t = Im(J '^-),t > is well defined |ItMK 65j. There is the 
well-known skew-product representation: 

log|Z t | + ^= / ^ = (B u + t lu ) ds , (25) 

where (B u + i , y u , u > 0) is another planar Brownian motion starting from log 1 -MO. Thus: 

H~ l = inf{t : H t > u} = / ds exp(2B s ) := A U {B). 

Jo 

For further study of the Bessel clock H, see e.g. |Yor80| . We also define the first hitting 
times T e c = inf{t : 6 t = c} and if 1 = inf{t :\9 t \ = c}. 



3.1 General results 

Bougerol's identity in law combined with the symmetry principle of Andre |And87j IGal08] 
yields the following identity in law (see e.g. |BeY12[ IBDY12a] ): for every fixed I > 0, 

H n = Ml) (26) 

where (77,/ > 0) stands for a stable (l/2)-subordinator. An example of this kind of 
identities in law is given for the planar Brownian motion case in the next Subsection. 
The main point in |BeY12| is that (126]) is not extended in the level of processes indexed 
by / > 0. 



3.2 Bougerol's identity in terms of planar Brownian motion 

Vakeroudis [VakllJ investigated Bougerol's identity in terms of planar Brownian motion 
and obtained some striking identities in law: 

Proposition 3.1. Let (/3„,m > 0) be a 1- dimensional Brownian motion independent of 
the planar Brownian motion (Z u , u > 0) starting from 1. Then, for any b > fixed, the 
following identities in law hold: 

i) H T i3 = T a(b) 11) e T p = c a{b) in) e T p = |c ( 6 )|, 

b y ' b b 

where Ca is a Cauchy variable with parameter A and 9 U = sup 0<s<u 6 s . 

Proof. 1) We identify the laws of the first hitting times of a fixed level b by the processes 
on each side of (JZ|) and we obtain: 

rpB (^) TJ 

I a(b) — n T^ 
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which is i). 

ii) It follows from i) since: 



e, 



(law) 



with (7 S , s > 0) a Brownian motion independent of (H u , u > 0) and (C u , u > 0) may be 
represented as {^ t b , u > 0) . 

in) follows from zi) again with the help of the symmetry principle. □ 

Using now these identities in law, we can apply William's "pinching" method [ Wil74[ 
McY82j and recover Spitzcr's celebrated asymptotic law which states that |Spi58|: 



2 (law) 
U t > Gi 



log t t- 



(27) 



with C\ denoting a standard Cauchy variable (for other proofs, see also e.g. |Wil74[ 
IDnr821 IMeY821 IBeW941 IYor971 IVaYllaj ). One can also find a characterization of the 
distribution of T c e and of tJ 6 *' in [Vakil] . First, applying Bougerol's identity ([ID in terms 
of planar Brownian motion we have: 



Proposition 3.2. For fixed c > 0, 



smh(6 c ) = P( T B) = 



(28) 



where N ~ A/"(0, 1) and the involved random variables are independent. 



Furthermore, we can obtain the following Gauss-Laplace transforms which are equiv- 
alent to Bougerol's identity exploited for planar Brownian motion: 



Proposition 3.3. For x > and m = 



cE 
cE 



7T 



2T e 



exp 



x 



7TT ( 



exp 



2T C * 

x 
2Tj e| 



Vl + X (C 2 + \0g\y/x + VT+X)) 
1 2 



(29) 



vTTx (y/l + X + x/i) m + {y/l + x - y/x) 1 



(30) 



Proof. For the proof of ( 129|) . it suffices to identify the densities of the two parts of (T58 
and to recall that the density of a Cauchy variable with parameter c equals: 



7r(c 2 + y 2 ) 

For (130|) . we apply Bougerol's identity with u = Tj 7 ' = inf{t : | | = c} and we obtain: 



sinh(S h) '*= j p (rr \o 



(law) (law) 
(Tj 91 ) = 



TP N . 



(31) 
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Once again we identify the densities of the two parts. For the left hand side, we use the 



following Laplace transform: for A > 0, E 



e 2 1 b 



i 



cosh(Afe) 



(see e.g. Proposition 3.7, 



p. 71 in Revuz and Yor |ReY99] ). We also use the well-known result [Lev80l IBiY87] : 

1 1 



E 



exp(i\B 



•hi. 



cosh(Ac) cosh(7rA|) 



-*(^' 1 ' dx. 



2tt cosh ( |) 



(32) 



Changing now the variables y = cx/tt, we obtain the density of B \ y \ which equals: 



(2coosh(f))-' 



, E ME ■ 

eye 2c + e 2c 



and finishes the proof. 



□ 



Vakeroudis and Yor in [VaYllaj IVaYllb] investigated further the law of these random 
times. 



3.3 The Ornstein-Uhlenbeck case 

Vakeroudis in [Vakil \ IVakl2j investigated also the case of Ornstein-Uhlenbeck processes. 
In particular, we consider now a complex valued Ornstein-Uhlenbeck (OU) process: 

Z t = z Q + Z t - A / Z s ds } (33) 
Jo 

where Z t is a complex valued Brownian motion, z G C (for simplicity and without loss 

of generality, we suppose again z = 1), A > and = T] e ' = inf {t > : \9f\ = c} 
(Of is the continuous winding process associated to Z) denoting the first hitting time of 
the symmetric conic boundary of angle c for Z. Then, we have the following: 

Proposition 3.4. Consider (Z*,t > 0) and (Uf,t > 0) two independent Ornstein- 
Uhlenbeck processes, the first one complex valued and the second one real valued, both 
starting from a point different from 0, and define T^ X \lI x ) = inf {t > : e xt Uf = 6} ; for 
any b>0. Then, an Ornstein-Uhlenbeck extension of identity in law ii) in Proposition 
\3.1\ is the following: 

where a(x) = argsinh(x) and Ca is a Cauchy variable with parameter A. 



Proof. First, for Ornstein-Uhlenbeck processes, is well known that |ReY99j . with (B t , t > 0) 
denoting a complex valued Brownian motion starting from 1, Dambis-Dubins-Schwarz 
Theorem yields: 

/ rt 
/ 1 i / „As, 



e- M I 1 + J e As dZ s 
e- xt (M at ), (35) 
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Let us consider a second Ornstein-Uhlenbeck process (U x ,t > 0) independent of the first 
one. Taking now equation (135|) for U x (1-dimensional case) we have: 



e M U t A = S (e 2M^, (36) 
where (5 t ,t > 0) is a real valued Brownian motion starting from 1. 

Second, applying Ito's formula to f )35|) and dividing by Z s , we obtain (at = jte 2Xs ds = 



2A r 

'dZ q \ , fdl 



Im j -7T- I Im 

hence: 

ef = el . (37) 

By inverting a t , it follows now that: 

T C W = ^ In (l + 2ATf *) . (38) 

Similarly, for the 1-dimensional case we have: 

Tl x) (U x ) = ±\n(l + 2\T*). (39) 

Equation fl3T|) for t = ^ In (l + 2AT b 5 ) , equivalently: a(t) = T 6 5 becomes: 
/iZ A _ /iZ A _ qB (^) ^ 

2f>(r/*) ~ jxH 1+2XT t) ~ 6 ~ a(fe) ' 
where the last equation in law follows precisely from statement ii) in Proposition 13.11 □ 



4 Multidimensional extensions of Bougerol's identity 

4.1 The law of the couple (sinh(^), sinh(Lf)) 

A first 2-dimensional extension of Bougerol's identity was obtained by Bertoin, Dufresne 
and Yor in |BDY12a| (for a first draft, see also |DuYll] ). With (L t ,t > 0) denoting the 
local time at of B, we have: 

Theorem 4.1. For fixed t, the 3 following 2-dimensional random variables are equal in 
law: 

(smh(B t ), sinh(L t )) ( '= } ((3 At ,exp(-B t ) X At ) {l = ] (exp(-fl t ) At ,X At ), (40) 
where (X u ,u > 0) is the local time of ft at 0. 
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Remark 4.2. Theorem \4-l\ can be equivalently stated as: for fixed t, the 3 following 
2- dimensional random variables are equal in law: 

(sinh(|£ t |),sinh(L t )) ( '= } (\(3\ At , ex V (-B t ) \ At ) {l = ] (exp(-ft) |/?k,A^). (41) 
Using now Paul Levy's celebrated identity in law (see e.g. \ReY99§ ): 

{{Bt - B t ,B t ),t > 0) (/ = } ((\B t \,L t ),t > 0) , (42) 
we can reformulate [JW or fcU\ )> an d we obtain: 

(sinh(fl t - B t ), sinh(5 t )) (/ = } (0 - 0) A „ exp(-S t ) (3 At ) 

{l = ] (exp(-fl t ) 0-P) At J At ). (43) 

The latter is particularly interesting when compared with the Wiener-Hopf factorization 
for Brownian motion. In particular, if we consider e q an independent exponential ran- 
dom variable of parameter q, then B Bq is independent of B Sq — B eq . This tells that the 
two random variables appearing on the right hand side of [Jcty , when taken at e q , are 
independent. 

Remark 4.3. Considering only the second processes of the first and the third part of fijffl 
(or equivalently of ffiW), we obtain a "local time" version of Bougerol's identity: 

sinh(L t ) { = } \ At , (44) 

which (as was shown in \BeYlty ). similar to the Brownian motion case, is true only for 
fixed t and not in the level of processes. 

Proof. ( Theorem \4-l\ ) 

From Remark 14.21 it suffices to prove ()41|) . 

First, we denote S p , p > an exponential variable with parameter p independent from B 
and g t = sup{w < t : B u = 0}. We know that (B u , u < g Sp ) and (^B gSp+u , u < S p - g Sp ^j 
are independent, hence Ls p and Bg p are also independent. We also know that L t and 
\B t \ have the same law. Hence, using the following computation: for every I > 0, with 
(ti,1 > 0) denoting the time L reaches I, 

P (L Sp >l)=P(S p > n ) = E [exp(-pn)} = exp(-/ y/ty), 

we deduce that the common density of L$ p and \Bg \ is: 

A/2pexp(— uy/2p), u > 0. 

Equivalently, we have: 

v^(|/3(l)|,A(l)) ( '= } (e,e'), 

where on the left hand side e and e' are two independent copies of Si independent from 

(3. 

For the second identity in law in Theorem 14.11 it suffices to remark that 
((3 At ,eM-Bt) A At ) {l = ] ( V / A/3i,exp(- J B t )v / A A x ), 
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and use a time reversal argument. 

For the first identity in law we use an exponential time S p and we compute the joint 
Mellin transforms in both sides in order to show that: 



2e(sinh(|i?| 5p ),sinh(L 5p )) ( '= } v / 2^(exp(- J B 5p )^ ^ p X At ). 

For further details we address the reader to [BDY12a| . □ 

Using now Tanaka's formula we can also obtain the following identity in law for 2- 
dimensional processes: 

Corollary 4.4. 

(45) 



(sinh(B t ),L t ) t > {l = ] (^exp(-B t ) (3 At , J exp(-B s )d\ Ai 



t>o 



where, in each part, the second process is the local time at level and time t of the first 
one. 

4.2 Another two-dimensional extension 

In this Subsection we will study the joint distribution of: 

(X«, X®) = (eM-Bu) £ d& exp(B„), exp(-2fl tt ) £ dtf> exp(2^)J , (46) 

where , v > 0), (£» , v > 0) and (B UJ u > 0) are three independent Brownian motions. 
Hence, we obtain a new 2-dimensional extension which states the following: 

Proposition 4.5. We consider (B^,t > 0) and {B^\t > 0) two real dependent Brow- 
nian motions, such that: 

d < B«, B® > v = tanh(5( 1) ) tanh(2£( 2 )) dv. (47) 

For the two-dimensional process (^Xu\Xu^, we have: 

(i) In the level of processes: 

(JffU?,«>0) { = ] (sinh(5«),isinh(25( 2 )), M >o) (48) 

(ii) For u fixed, 

( X W X®) ( B (l) B (2) \ (49) 
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Proof. Let us define: 

X^=exp(-aB u ) d^ exp(aB v ), (50) 
Jo 



where a = 1, 2. By Ito's formula, we have: 



*L a) = + J^ex V (-aB v )(-adB v ) + ^ex P (-aB v )dv^j (jT d^ ex P (aB h )^j 



Hence: 



*?> = & - dB„ X? + \ £ XV to 

- £*pJ{ i+ W) +l 2£ x ? ) *'- (51) 

and 

Xi 2) = e - 2 / dB v + 2 / dv 
Jo Jo 

= [ dvl 2) \j (l + 4 (XP) + 2 jf X <*> d„, (52) 

where (r^\v > 0) and (rj^\v > 0) are two dependent Brownian motions, with quadratic 
variation: 

d < r/ 1 ), ^ >„= 2Xv X v . dV (53) 



1+ (X«)^( 1 + 4(XP) 2 ) 
Thus, we deduce that the infinitesimal generator of is: 



«9 2 . QN d 2 d 2 



> + x2 i)t^ + ( 1+ Ax D 7^2 + 



x\ d d 



+ 2x 2 —. (54) 



2 dxi dx 



dx\ ' dx\ dx\dx2 

Let us now study the couple: 

(a^f) = (sinh(5; i) ),isinh(25; 2) )^) , (55) 

where (Bf~\t > 0) and (Bf\t > 0) are two dependent Brownian motions. By Ito's 
formula we have: 

x? = sinh(5«) 



[ cosh(B«) dB^ + \ f sinh(flW) dv 
Jo ^ J 

l + (xP)A dB^ + \ tx^dv, 



(56) 



'o 
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and: 



xf ] = ^smh^f) 



cosh(2^ 2) ) dBi 2) + / sinh(2^ 2) ) dv 

JO 



= J* y (l + 4(xl 2) ) 2 ) dB^ + 2 jf x< 2) dv. (57) 

Moreover, using ( |47|) : 

d < Binh(sW), - sinh(25^) = cosh(BW) cosh(2£( 2 )) d < B^\ B^ >v 

= 2sinh(£( 1) ) ^sinh(2^ 2) ) dv. (58) 

Finally, we have that (xu\ x^) has the same infinitesimal generator with (^Xu\Xu^. 
Hence, we get part (i) of the Proposition. 

For part (ii), we fix u and we have: 

( f sii 1 h(Bi 1 )),-Bmh(2S«)V = ) (fil , «,P?L , A (59) 

^ V U h 2 V " 7 \ U exp(2B„))'^(/ u ^ exp(4S„))y ' V ; 

where , u > 0) and , f > 0) are two dependent Brownian motions and (B v , v > 0) 
is another Brownian motion independent from them. Now, from (l59p . we obtain (IBDj) . □ 

Remark 4.6. From / f5Pj) . toi£/i p u (x,y) denoting now the density function of the couple 
(sinh( J Bi 1) ), isinh(2M 2) )), we have: 



Pu(x,y) = E 



1 exp (-x 2 /2 £dv exp(2B v )) exp (-y 2 /2 £ dv exp(4B„)) 
2^ 



Jq dv exp(2B v ) J J Q " dv exp(AB v 



(60) 



In theory, we should be able to compute this probability density as we know the joint 
distribution of the couple of exponential functionals (see e.g. fAMShOlf ). 

4.3 A three-dimensional extension 

Alili, Dufresne and Yor, in |ADY97| . obtained a 3-dimensional extension of Bougerol's 
identity: 

Proposition 4.7. The two following processes have the same law: 

e Bt [ e- B "dp u ,B t ,p t ;t>o\ {l = ] {smh(B t ), B' tJ G' t ;t > 0} , (61) 



o 

where: 

' .III- r<< r> ^ ( 62 ) 



Gi = fod§i^-f>MB s )dG s 



with (G t ,t > 0) denoting another Brownian motion, independent from B. 
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Remark 4.8. We remark that with: 



tanh(x) 



we have: 



and 



is a 2- dimensional Brownian motion. 



Proof. ( Proposition \4- 7\ ) 

First proof: Using Ito's formula, we deduce easily that each of these triplets is a Markov 
process with infinitesimal generator (in C 2 ' 



1 ^ 2 d 2 + 1 d 2 + 1 d 2 + ^ d 2 + d 2 + ^ d 

2 dx 2 2 dy 2 2 dz 2 dxdy dxdz dx 

The proof finishes by the uniqueness (in law) of the solutions of the corresponding mar- 
tingale problem. 

Second proof: First, we admit that the identity in law is true. Then, if we replace 
on the left hand side (B s ) by (B' s ) and (j3 s ) by (G' s ), we have necessarily: 



sinh(^) (/ = } e^ te^dG't, 
Jo 



(66) 



which is essentially a (partial) inversion formula of the transformation ( l64jh 

Equation (166]) can be proved by using Ito's formula on the right hand side. □ 

Gruet in |ADY97j also remarked that: 

Proposition 4.9. There exist two independent linear Brownian motions V and W and 
a diffusion J starting from satisfying the following equation 

dJ t = dW t + ^ tanh( J t )dt, (67) 



such that, 



Hence, the two following 3-dimensional processes: 

(«cp (B t + jT exp (-B B - |) d(3 s , B u ft; t > 

and 

(smh(J t ),B t ,p t ;t>0), 

are equal. 
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Proof. This result follows from a geometric proof and it is essentially an explanation of 
the second proof of Proposition 14. 7\ at least for v = 0. For this purpose, we can compare 
the writing of a hyperbolic Brownian motion in the half-plane of Poincare, decomposed 
in rectangular coordinates with the equidistant coordinates |Vin93] . For further details, 
see the Appendix in |ADY97| due to Gruet. □ 



5 The diffusion version of Bougerol's identity 
5.1 Bougerol's diffusion 

Bertoin, Dufresne and Yor in a recent work |BD Y12b| generalized Bougerol's identity in 
terms of diffusions. First, we remark that from Proposition 12.11 we have that (see also 
(ADY97J): 

(smh(B t ),t > 0) (/ = } (exv(-B t )/3 4 o h t > o) . (69) 

In particular, using Lamperti's relation (see e.g. |Lam72| or |ReY99| ) we can invoke a 
Bessel process R^ s ' independent from B in order to replace the right hand side of (169|) by: 

exp {-B^)R^,t>0 

which turns out to be a diffusion (named Bougerol's diffusion) with a certain infinitesimal 
generator. Hence, we obtain the following: 

Theorem 5.1. With Z = Z^ and Z' = Z^ 5 '^ denoting two independent squared Bessel 
processes of dimension 8 = 2(1 + //) and 8' = 2(1 + v) respectively, starting from z and 
z' , the process: 



X t ee ee exp (-2flW) Z A ^ = ^ 



t > 0, (70) 



u=A<l v) 



is a diffusion with infinitesimal generator: 

2x(l + x)D 2 + (8 + (4 - 8')x) D, (71) 

where 8' = 2(1 + v). 

Remark 5.2. There is a discussion in jJaWlty concerning the particular case where 
the diffusion with generator given in IF/1\) is the hyperbolic sine of the radial part of 
a hyperbolic Brownian motion ( or equivalently the hyperbolic sine of a hyperbolic Bessel 
process) of indexa G (—1/2, oo) (see fJaWHty Theorem 2.25, formula (46), p. 15). In that 
case, with R t denoting this hyperbolic Bessel process starting from x andY t = e Bt ~( a+1 / 2 ^ , 
for any w > 0, t > 0, 

(sinh(^), t > 0) ( '= ) (K^Xy**, t > 0) , (72) 
where S is a Bessel process of dimension 2(1 + a) independent of B, and Sq = smh(x). 
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Proof. ( Proposition 15. ij) 

Applying Ito's formula to the process X, we obtain: 

X t = J\xp (-2SM) d (Z^) + J Z^d (exp (-2£#>)) . (73) 

For the second integral in (|73|) . Ito's formula once more yields: 

d (exp (-2BM)) = -2 exp (-25^) + i/du) + 2 exp (-25^) 
= -2 J X u dB u + 2iy-u^ J X u du . 

° AS 1 ° 

Thus: 

/ Z A ^d(exp(-2B^)) = -2 J X u dB u + 2(J_-j^J X u du . (74) 

° ° 4-5' ° 

For the first integral in (173|) . we recall that, with 7 denoting another Brownian motion 
independent from B (thus independent also from Z): 

dZ s = 2y/z~ s d% + 5 ds . (75) 

Hence: 



dZ A ( V ) = 2JZ A (v)dj A (v) + 5 exp (2B^) du 



2JZ aM exp (BM) d<y u + ^ exp (2B^) du , (76) 



with 7 denoting another Brownian motion, depending on 7 and on B. 

The proof finishes by some elementary computations from (173 p . using (174)) and (TT6"|) . 

Finally, using Lamperti's relation, which states that: 



exp (2fl t w ) = Z' A ^ , (77) 



we obtain the last identity in fPTO]) . □ 

We may continue a little further in order to obtain the following result relating the 
diffusion X with its reciprocal (recall that: = Jq ds exp 

Corollary 5.3. The following relation holds: 

' X™ . (78) 



x^ rojp 

-™-H. 



Proof. It follows easily by some relations involving the changes of time: 

A[ u) = [ ds exp (2Bi u) ) ; = [ ds exp (2B^) ; (79) 

Jo Jo 
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Moreover, we remark that [H[ '.) is the inverse of (A^) and (H^) is the inverse of (A^). 
We also need to use: 

H {u,n) = j U ds . H M = f U ds 



Simple calculations now yield: 



o xr H ' Jo x<f' v) ' 



= H% , (82) 



Finally, using (T70|) we have: 

H % = I ' (83) 

and we obtain easily the result. □ 



5.2 Relations involving Jacobi processes 

In this Subsection, we present a particular study of Theorem 15.11 in terms of the Jacobi 
processes Y ' as introduced in Warren and Yor |WaY98] (see also the references therein 
for Jacobi processes), due to Bertoin, Dufresne and Yor [BDY12bJ. First, we recall some 
results involving Jacobi processes: 

Proposition 5.4. (Warren and Yor |WaY98] . Proposition 8) 

With T = inf{u : Z u + Z' u = 0}, there exists a diffusion process (Y u = Y^' s , u > 0) on 
[0, 1], independent from Z + Z' such that: 

-J^ = Y U<T . (84) 

Z u + Z' u J o z^zZ 

We remark that Y' = 1 — Y is the Jacobi process with dimensions {S',5), and Y has 
infinitesimal generator: 

2y(l-y)D 2 + (5-(5 + 5>)y)D . (85) 
Now, X defined in (17U1) and Y can be related as following: 
Proposition 5.5. The following relation holds: 

Y w 

= X rw dv = X rw dv , (86) 

1 - Y„ j Y> JO 1-Y V 



or equivalently: 

Y, 



V "W 

At 



1-Y W 



17) 



k dv 



1+X V 
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Proof. ( Proposition \5.5\) 
First, from fPTUl). we have: 



Conversely, 



where Hu^ = J Q " is the inverse of A^ u \ However, using the Jacobi process Y, 



u=A 



Z u 

~zj. 



11 II 



Y 

■*- in 



Z' u 



Zs+Z' s 



and moreover: 



ds 



nu dv 
Yf, 



(89) 



(90) 



/„ (Z s + Z' s ) (l-Y Hs ) 

Plugging now f[9"0"j) to ( THH|) and comparing to f lS^j) . we obtain (jSSj) . For (I8"T|) . it suffices to 
remark that k — > ls the inverse of the increasing process w — » J™ ^p- . □ 



6 Bougerol's identity and peacocks 

Hirsch, Profeta, Roynette and Yor in [HPRYllJ, studied the processes which are in- 
creasing in the convex order, named peacocks (coming from the French term: Processus 
Croissant pour l'Ordre Convexe, which yields the acronym PCOC). Let us first introduce 
a notation: for W and V two real- valued random variables, W is said to be dominated by 
V for the convex order if, for every convex function ip : IR — > R such that £7 [ [ ^ (T^) | ] < 00 
and i£[|^(V)|] < 00, we have: 

Ety(W)\ < Ety{V)l (91) 

(c) 

and we write: W < V. 

(c) 

A process (G t ,t > 0) is a peacock if, for every s < t, G s < G t - Kelleler's Theorem 
now (see e.g. |Kel72|. IHPRYllj IHiR12] ) states that, to every peacock, we can associate a 
martingale (defined possibly on another probability space than G). In other words, there 
exists a martingale (M t ,t > 0) such that, for every fixed t > 0, 

G t ( = ) M t . (92) 

The main subject of [HPRYllJ is to give several examples of peacocks and the associated 
martingales. 

We return now to Bougerol's identity and we remark that (see also |HPRYli~] . paragraph 
7.5.4, p. 322), for every A > 0, (smh(XB t ), t > 0) is a peacock with associated martingale 

(A j*e x ?°d ls ,t >o) (see e.g. ©). 

Moreover, for every A real, ^e~V* sinh(A-B t ), > 0^ is obviously a peacock, as it is a mar- 
tingale. This is generalized in the following: 
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Proposition 6.1. ( |HPRYllj . Proposition 7.2) 

The process (e M * smh(\B t ), > 0) is a peacock if and only if fi > —\- 



Proof, i) First, we suppose \± > — 4^-. Then, for s < t: 

/ . . A 2 A ( c ) ,.. , A 2 x / . . A 2 



e^sinh(AS t ) = e {/1+ V)M smli^)^*) > e^+VW sinh(A5 t )e-^*J 
> e (/1+ ^ )s (sinh(A5 s )e-^ s ) . 
zz) Conversely, Ito-Tanaka's formula yields: 

£ [| sinh(AB t )|] = E [sinh(A|B t |)] = e^*A / 



/o V2 



hence: 



S[|e^sinh(//S t )|] Ae^+^A / 



e 2 



v/2 



7T.S' 



which means that if [i < x 



2 ' 

E [|e^sinh(/i^)|] 'A ro . 

However, x — > \x\ is convex and if (e Mi sinh(/iS t ), t > 0) was a peacock, then E [|e M * sinh(//S t ) 
would increase on t, which is a contradiction. □ 



7 Further extensions and open questions 

In this Section, we propose some possible directions to continue studying and possibly 
extending Bougerol's celebrated identity in law (for fixed time or as a process). 

First, the natural question posed is wether this identity can be extended to higher 
dimensions. This very challenging question has already been attempted to be dealt with, 
and in this paper we've presented several extensions, at least for the 2-dimensional (and 
partly for the 3-dimensional) case. 

Another natural question is wether we can generalize Bougerol's identity to other 
processes. For this purpose, we may think in terms of a diffusion, as introduced in 
Section [5j It seems more intelligent to start from the right hand side of ([1]) and try to 
see, e.g. in (170 p . for every particular ratio of processes, which is the corresponding process 
on the left hand side (this process could be named "Bougerol's process"). 

In particular, it seems interesting to investigate a possible extension in the case of Levy 
or stable processes. To that end, we could replace the ratio of the two squared independent 
Bessel processes in f lTDj) by e.g. the ratio of two exponentials of Levy processes, and 
investigate the process obtained after the time-change. However, this perspective is not 
in the aims of the present work. 

Finally, another aspect which could be further studied is the applications that one 
may obtain by the subordination method, as presented in Section [3j Following the lines 
of this Section, one may retrieve further results and applications, others than for the 
planar Brownian motion case (see also |BeY12t lBDY12a"] ). 
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A Appendix: Tables of Bougerol's Identity and other 
equivalent expressions 



Using now the notations introduced in the whole text, we can summarize all the results 
in the following tables (u > 0, wherever used is considered as fixed). 



A.l Table: Bougerol's Identity in law and equivalent expressions 
(u > fixed) 

With a(x) = argsinh(x), and B, (3 denoting two independent real Brownian motions, for 
u > fixed, we have: 



(law) 



sinh(£ u ) P {A 

u {B)=£dsex P {2B s )) (Bougerol's Identity) 



2) sinhd^l) {l = ] \P\ (Au{B)) 



(law) 



3) smh(B u ) = P( Au (B)), B u = sup s < u f3 s 



4) E 



^2wA u (B) 



A. 2 Table: Bougerol's Identity for other 1-dimensional processes 
(u > fixed) 

We use fi, v reals and we define: 

B^ = B t + f it, f3^ = (3 t + ut, 



Af ] = fds exp(2SM), 
./o 



e: a Bernoulli variable in {—1, 1}, (i?t,£ > 0) a 2-dimensional Bessel process started at 
0, S an arcsine variable, and > 0) a diffusion with infinitesimal generator: 

2^ + ^ tan%) + co^%)J^' 

starting from y = argsinh(x). _B^, e, S and i? are independent. Then, for w > 
fixed: 
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5) (smh(Y t M ),t > o) ( =" ) (exp^) (x + ^exp(-B^)d^ ,t > o), a: fixed 



6) sinh(Yj^) {l = ] £exp(B^)d^\ 



7) sinh(_B M + e t) = J"exp(_B s + s)d(3 s , e: Bernoulli variable in { — 1, 1} 



8) /3 a m (l = ) (2S - 1)0 (flM, /r 2 + (#) 2 ) , = ^cosh^-e 2 *-!, z > \x 



A. 3 Table: Bougerol's Identity in terms of planar Brownian mo- 
tion (u > fixed) 

We define (Z t , t > 0) a planar Brownian motion starting from 1. Then t = Im( f Q * ^§ s -),t > 
is well defined. We further define the Bessel clock H t = J* * jjjp = A~ l {B) and the first 

hitting times: T e c = in£{t : 9 t = c} and tJ 9 ' = inf{t : \0 t \ = c}. Then, with (C c , c > 0) a 
standard Cauchy process, a Cauchy variable with parameter y, N ~ jV(0, 1), (Z t A , t > 
0) and (Uf,t > 0) two independent Ornstein-Uhlenbeck processes, the first one complex 
valued and the second one real valued, both starting from a point different from and 
Tl X) (U x ) = inf {t > : e xt U t x = b}, for 6, c> fixed: 
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9) sinh(C c ) ( '= } /3 (T *) ^ y/TfN, c > fixed 



(law) 



10) H T p=' T*L, Tf = mi{t : B t = y}, a(x) = argsinh(a;), b > fixed 



(law) 



11) 9 T fj = C a (b) , b > fixed 



12) 9 T ft ( = ) |C a ( 6 )|, 9 U = sup s < u # s 



13) E 



exp l — 



l+x 7r(c 2 +log 2 ( v / z+ v / T+i)) 



, b > fixed, x > 



14) £ 



(s Us) 



( %/l+X+ %/x) < + ( VT+x' - y/x) < 



,x>0,( = f c 



15) 9 Z T l ){uX) {l = ] C a(b) (OU version) 



A. 4 Table: Mult i- dimensional extensions of Bougerol's Identity 

In the following table, (L t ,t > 0) and (X t ,t > 0) denote the local times at of B, j3 
respectively and: 

WUP) = (eM-Bu) £de^MB v ), exp(-2B u ) £ df® exp(2B v )) , 

where ($\v > 0), ($\v > 0) and (B u , u > 0) are three independent Brownian motions. 
Moreover, we denote by (B^\B^) and {(5^,(5^) two couples of dependent Brownian 
motions (independent from B), such that: 



d < BW,BW > v = tanh(5«) tanh(25^) dv, 



and, for w > fixed: 
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Finally, 

with (G t , t > 0) denoting another Brownian motion, independent from B and J a diffusion 
starting from satisfying: d.J t = dW t + \ tanh(J t )dt, where W stands for an independent 
Brownian motion. Hence, for u > fixed: 



16) 


, . (law) 

sinh(L u ) = \ Au 


17) 


(sinh(5 u ),sinh(L u )) (, = } (f3 Au , exp(-5 u ) X Au ) ( '= } (exp(-£ u ) f3 Au ,\ Au ) 


18) 


[smh(\B u \) sinh(L u )) ( '= } exp(-£„) X Au ) (l = ] (exp(-S u ) |/?k,A A J 


19) 


(sinh(5„ - J B tt ),sinh( J B u )) ( '= } ((£ _ /3) Au , exp(- J B u ) {l = ] (exp(- J B u ) (0 - /3)a u , AO 


20) 


(lain} / / \ 

(sinh(S t ),L t , i > 0) [ = ) (eM-Bt) f3 At , J* exp(-B s )d\ As , t > Oj 


21) 


(* t (1 U t (2) ,i > Oj ( = W) (sinh( J B t (1) ),isinh(25f ) ),t > Oj 


22) 


y(2)\ (iotu) / (1) (2) \ 

v ' u y ~ l v p (/ ( r^cxp(2B„))' /3 (/ l r^cxp(4^ ) ) y ) 


23) 


(e B < £ e^/3 tt , 5,, A; t > 0) ( '= } (sinh(S t ), B' t , G> t ; t > 0) 


24) 


(exp (B t + |) £ exp - f ) d(3 s , B t , fa t > o) ( l = ] (sinh( J t ), 5 t , A; t > 0) 
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A. 5 Table: Diffusion version of Bougerol's Identity (relations in- 
volving the Jacobi process) 

Let Z = and Z' = Z^ be two independent squared Bessel process of dimension 
S = 2(1 + fi) and 8' = 2(1 + v) respectively, starting from z and z', and X t = X^'^ a 
diffusion (named "Bougerol's diffusion"), with infinitesimal generator: 

2x(l + x)D 2 + (5 + (4 - 8')x) D, 

and Y = Y s > 6 ' the Jacobi process. Then, for t,w,k > 0: 

25) ^r^exp(-2 B «)z 4 ,„, = |j|^ M 



26) -. ^ V- — X rw dv_ — X rw dv 

L-y w j y ^ j 1 _ Yv 



27) X k = ^ 
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